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On the generation of sequential unitary gates from continuous time 
Schrodinger equations driven by external fields 
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In all the various proposals for quantum computers, a common feature is that the quantum circuits 
are expected to be made of cascades of unitary transformations acting on the quantum states. A 
framework is proposed to express these elementary quantum gates directly in terms of the control 
inputs entering into the continuous time forced Schrodinger equation. 
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I. INTRODUCTION 



In quantum information, the "computing" with quantum states is accomplished by applying sequences 
of discrete unitary gates, i.e. predetermined elements of the transformation group that acts on the state 
^ , of the quantum system and determines its dynamical properties. Such set of unitary transformations is 

' normally specified in terms of the corresponding infinitesimal generators and, in order to guarantee arbitrary 

manipulation of the quantum states, it should be enough rich to generate the whole transformation group. 
In quantum computing, this property is called universality of the gates and in j|] it is shown that "almost 
^— ^ ' all gates are universal" , see also |2y, for a control theoretic reformulation of the same idea. The problem 

, is treated also in Q| by recollecting results on the uniform generation of the rotation groups. The main result 

' is the following: arbitrary state transfer is achievable by means of a number of unitary gates equal to the 

dimension of the transformation group (henceforth SU{N) of dimension N'^ — !)■ Such unitary gates do not 
^ • have to correspond to linearly independent infinitesimal generators: in fact, owing to the semisimple nature 

" ' of SU{N), almost all pairs of infinitesimal generators in the Lie algebra su(A^) allows to achieve arbitrary 

■ state transfer by cascading N'^ ~ 1 exponentials of the two generators in alternate sequence. For the single 
' qbit, this is the idea behind for example the ZYZ parameterization of rotations of the corresponding group 

^ . of transformations SU{2). For TV qbit systems, parameterizations aiming at the same scope are proposed 

for instance in [T^. From a technological viewpoint, this is a great simplification as it means that in 
practice arbitrary manipulation can be done by sequencing control pulses applied along (almost) any pair of 
laboratory fixed axes. 2D NMR spectroscopy is one field where such techniques of sequencing pulses from 
. , two different directions has reached the level of a very sophisticated science ^ . The cascade of unitary 

■ transformations obtained in this way is a product of exponentials, each exponential containing a parameter 
expressing the time duration of the pulse (scaled by a constant depending on the "strength" of the applied 
field). By varying these N'^ — 1 parameters, every quantum state is reached. 

In general, the philosophy of control synthesis for states that evolve in continuous time governed by 
differential equations (like the Schrodinger equation) is to look for input-parameterized solutions of the 
dynamics and then to select the control inputs in order to accomplish the given tasks (for example state 
transfer). As mentioned above, in the sequential model for quantum computing it is normally expected that 
the elementary gates that constitute the "quantum circuits" will be made of discrete sets of unitary operators 
corresponding to fixed values of the N"^ — 1 parameters mentioned above, applied in a fixed order. In other 
words, what in a different context would be simply a parameterization of a group manifold, like a set of Euler 
angles, in the context of quantum computing becomes the "hardware" basis for the construction of quantum 
circuits. This "mimicking" the behavior of classical circuits is motivated by the need to simplify as much as 
possible the influence of the dynamics in the logic of the circuits and of the "hardware" requirements of the 
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circuits, but it is complicated by the continuum of values in which the quantum state can exist, as opposed 
to its classical counterpart (i.e. the bit). 

The main scope of this paper is to provide a framework for reconciliating the continuous time forced 
Schrodinger equation with the "discrete dynamics" of quantum computers. The formalism used is that of 
the Wei-Norman formula, which relates the Magnus expansion (suitable to represent the continuous time 
evolution of the Schrodinger equation) with the product of exponentials expansion corresponding to complete 
sets of elementary gates. Essentially, both expansions define local coordinate charts on the group of unitary 
transformations by means of a basis of the corresponding Lie algebra and via the exponential map. That 
is why these are sometimes referred to as exponential coordinates (or, in some cases, canonical coordinates). 
For compact groups, the two corresponding parameter spaces are different: one is a solid sphere and the 
other a cube, and the two sets of coordinates are related by nonlinear differential equations. An algorithm 
for the explicit computation of such differential equations (i.e. of the Wei-Norman formula) in terms of the 
structure constant of the Lie algebra was recently proposed in |^] . As for all series expansions on semisimple 
Lie groups, existence is not a global property, see ||2^. In fact also the Wei- Norman formula has a singular 
locus which is an algebraic set of the parameter space. 

Another control-theoretically sound way to design inputs that resemble elementary quantum gates is to 



use piecewise constant controls |16| . The methods proposed here could be seen as the differential version 
of the control design obtained via piecewise constant controls. In fact, also the flow induced by piecewise 
constant controls looks like a product of exponentials which resembles a cascade of elementary gates. 

For sake of comprehension, a parallel with what happens in the control of robotic chains |^ is useful to 
understand this relation: while at each time instant one can express the position of the end-effector of a 
robot in terms of joint angles and topology of the robot via a static map, called forward kinematics, (in 
quantum control: a map constructed from the Lie group decompositions of the unitary propagator, i.e. from 
products of exponentials), in order to describe dynamic changes it is convenient to pass to the differential 
forward kinematics i.e. to the Jacobian map obtained by differentiating the static map mentioned above. 
In Robotics, the use of differential maps is crucial whenever a closed form expression of the inverse of the 
forward kinematics (giving the values of the joint angles in terms of the end effector position) is not available. 
The Wei-Norman formula is exactly the Jacobian of the change of coordinates between single exponential 
and product of exponentials. Furthermore, as we will see, it automatically allows to bypass what is the main 
problem of the piecewise constant control methods, that is to say finding an explicit value of the — 1 
parameters that constitute a decomposition of an element in SU{N) for > 2 (in the robotic case this would 
correspond to the inversion of the forward kinematic map). Lastly, while the Lie group decompositions used 
in piecewise constant control methods |l^ are not directly related to the positioning of the laboratory 
equipement that produces the control fields on the quantum system, this factor is implicitely considered in 
the product of exponentials we will use here. 



II. EXPONENTIAL COORDINATES ON SU{N) 

Consider a closed finite level quantum system described by a state 1-0) evolving according to the time 
dependent Schrodinger equation 

th\m) = H{tm = (i/o + H,{t)) m)) 
\m) - i^o) 

where the state of the A^-level quantum system 1-0) lives on the sphere in A^-dimensional complex Hilbert space 
§^~^ = {10} G s.t. ('0|'0) = ijj and the traceless Hermitian matrices iJp and Hi{t) are respectively 
the constant internal Hamiltonian and the external time-varying Hamiltonian, this last representing the 
interaction of the system with the control fields. The solution of (|l|) is normally written in terms of the 
unitary propagator U{t) £ SU{N): 
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with U{t) satisfying an equation similar to (||) but lifted from the sphere ^ to the special unitary group 
SU{N): 



inU{t) = H{t)U(t) = {Ha + Hi{t)) U(t) 
C7(0) = / 



(2) 



From now on, we shall use atomic units h = \. Assume Ai, . . . , An, n — iV^ — 1, are skew-hermitian matrices 
forming a basis of the Lie algebra s\x{N) and that the Hamiltonian H{t) can be written in this basis as 



-iH{t) = -I {Ho + H,)^J2 "^3^0 + 

where Oj and Uj = Uj{t) are respectively the (constant) components of the free Hamiltonian along the basis 
directions and the (time- varying) control parameters of the interaction part, some of which might be zero if 
the control along the corresponding direction is missing. 

If T is the Dyson time ordering operator, the solution of (^ is formally written in terms of the exponential 
map as 



U{t) = Texp (y-i j^^ H{T)dT^ = Texp ^-i J Y^iaj + Uj{t))Aj 



(3) 



normally referred to as Magnus expansion |^ . Alternatively, instead of (^) , the solution of (||) locally admits 
an expression in terms of product of exponentials, first due to Wei-Norman 



U{t) = exp(7iyli) . ..exp{jnAn) 



(4) 



with generally time dependent parameters ji — ^i{t). The relation between the two expansions ^ and (^) 
is given by the so-called Wei-Norman formula, which expresses the functions ji {t) in terms of the -I- Ui {t) 
via a system of differential equations: 



S(7i, . . . ,7„) 



71 




fll + Ul 









7.(0) = 



(5) 



with the n X n matrix 2 analytic in the variables 7^. Since the adjoint maps can all be written as (summation 
convention over repeated greek indexes): 



3) Ml A'2 



rf^' A 



(6) 



/=0 (=0 

the matrix S of elements (S)fei = is defined in terms of the ji and of the structure constants as: 



3 = 1 



m — 1, 



,n 



(7) 



When the cardinality of the basis ordering is followed in (^, the matrix S assumes also the meaning of map 
between canonical coordinates of the first kind and canonical coordinates of the second kind, see |2l| . In this 
case, since 7i(0) = 0, S(0) = / and thus S is locally invertible. However, any fundamental parameterization 
of SU{N) could be used: various possible sets of Euler-like angles for SU{N) are discussed in |l^, 
For example, we will concentrate on the ZYZ parameterization of SU{2) in the example of Section ^V| . 
Because of the semisimplicity of SU{N), all parameterizations lead to a Wei-Norman formula that is subject 
to singularities and as such has only a local validity. Call S the singular set of S. By inverting S when 
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possible, equation ^ assumes the more traditional aspect of a system of first order differential equations in 
the 7i variables: 



71 



In 



fll + Ul 



7^(0) =0 



(8) 



If the time evolution of one of the two vectors of coordinates "fi {t) or + Ui (t) is known, the formulas (^) or 
can be used to obtain the other one. A method for the explicit closed-form calculation of the Wei-Norman 
formula is proposed in |^ and will not be repeated here. 



A. Analysis of the parameter space of the two expansions 

While (I) is global, (|) is valid only as long as det(^) ^ 0, and thus the nonsingularity of ^ needs to 
be checked at the point of application. In general the exponential map does not posses any good property 
globally: for example for semisimple Lie groups it is surjective but not injective. The Lie algebra 5u(A'^) being 
"compact" implies that the parameter space F C M" of the coordinates corresponding to the exponential 
map, i.e. the set of values of the real coefficients needed to cover the whole Lie group under the exponential 
map, is a bounded domain, corresponding to the principal values of the logarithm map. For the single 
exponential representation and the basis Ai chosen above, the parameter space is a solid sphere of radius 
27r. 2tt is a conventional choice, see 0, p. 127, which originates from the standard choice of a factor i in 



the exponentiation of the Pauli matrices, convention that we also follow in the example of Section IV. The 
radius could obviously be rescaled by any real normalization factor. In fact, given any vector A e su(A^), the 
straight line through the origin 7^, 7 G M, is mapped to a one-parameter subgroup (a geodesic curve if the 
metric is obtained from the Killing form) of SU(N). Because of compactness, exp is periodic and for such 
straight line it holds that — gi^ for some real constant r depending on the norm of A. Since this 

is true for all A = '^j^j^ the principal values of the logarithm map will correspond to a solid ellipsoid 

in general, and to a solid sphere in M" of radius 2tt when the Aj are suitably normalized and the aj are such 
that X)j=i = 1- 

The lack of global properties of the exponential map is amplified in the products of exponentials such as 
(^. Call pexp the homomorphism 

pexp : F SU{N) 



7 ^ n ^ 

k=l 



where p{k) represents the ordering of the basis elements corresponding to the parameterization chosen for 
SU{N) (in @ p{k) = k, in the ZYZ parameterization of 31/(2) used below p(l, 2, 3) = 3, 2, 3). Although 
we have 7fc G M, also in this case the parameter space, i.e. the subset of ffi." containing the principal values 
of pexp^"'^, is a bounded domain in R. However, since each 7^ enters into pexp(7) independently from 7^-, 
k j, this time the parameter space Fp C M" corresponds to a cube in R" centered in the origin and with 
sides of length An for the same basis Aj for which the logarithm has domain of principal values equal to the 
sphere of radius 27r. 



III. ELEMENTARY QUANTUM GATES AS OUTPUT OF A SYSTEM OF DIFFERENTIAL 

EQUATIONS 

Calling 7 = [71 . . . 7n]^, a = [ai . . . an]"^ and u = [ui . . . Un]'^ , from (|[) and (^ the state of the quantum 
system can be thought of as the output of the following system of first order nonlinear differential equations: 

7 -s-i(7)(« + ^) 7(0) = 7erp\E ^^^^ 

\ip) = pexp(7)|-0o) 
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While the differential equation (^) was linear in U{t) and affine in the controls, the differential equation ( [lO| ) 
is still afhne in control but it is highly nonlinear in the variables 7. 

The study of systems of nonlinear differential equations like (10) is quite common in systems theory, 
provided we call 7 the "state vector" and the "output". As they are quite confusing in the present 
context, these (and other) systems theory notations will always be recalled between quotes. 



A. Coherent states and 7-coherent states 

In the whole paper we will always neglect global phase factors e"^ as they are irrelevant for our purposes: 
|V') = e^1V')- 

In general, an action of a Lie group on a manifold is said transitive if each pair of elements of the 
manifold can be joined by an element of the Lie group. On the Lie group SU{N) acts transitively: 

V Itpo) G there exists a unitary transformation U £ SU{N) such that = C/|?/'o)- The manifold 

in this case is a homogeneous space of the Lie group. 

The following definitions are standard for a quantum mechanical system (see [ p^ p^). The maximal 
isotropy subgroup of \ipo) is given by 

- {/i e SU{N) s.t. h\iPo) = Ho)} 

Since the action is transitive, H^^^^ is a subgroup. In fact, if hi, h2 G ^|j/jo)' then ^2^i|'0o) = ^2|'0o) — l^o) 
and h^^lipo) — \^o). States that differ by an element are indistinguishable from |^o): if ?7 G H^.^i,^^ 

then \ip) = Ulipo) = \^o): again (and always) up to an irrelevant global phase factor. Calling C|0p) the coset 
space of IV'o) in SU{N), i.e. the space of equivalence classes determined by the isotropy subgroup, 

each U G SU{N) can be decomposed into U = ilh with G q^o) and h G H\^^y The "effective" change in 
^/'o) given by a unitary transformation U = Qh corresponds to the coset component alone: 

W ^Ul^Po) =nh\i;o) =n\i;o) 

The set 

{r!|^o) s.t. r!Gq^„)} -q^„)|v^o) (12) 

is called the set of coherent states of j-^o)- As can be seen from (^2|), the coherent states of j-^o) are obviously 
in 1-1 correspondence with the elements of the coset space q^,,). 

For the product of exponentials map (^, the preimage pexp~^(?7) may vary from point to point. In 
particular, pexp~"'^(/) varies according to the order p{k) selected. For example, in the canonical coordinates 
of the second kind {p{k) — k) pexp^^(/) = {0}, but in general it might be a nontrivial set. Call such set 
Hj{U), the 7-isotropy subgroup at U of the map pexp at [/ G SU{N) 

H^{U) = {7 e Fp s.t. pexp(7) = U}= pexp-i(C/) 

For example, H^{I) is the set of all coordinates 7 G Fp that certainly do not produce any effect on 
because their corresponding unitary operator is the identity. More generally, all 7 such that pexp(7) G ^^|i/;o) 
will produce no effect on an initial quantum state IV'o)- We indicate the composition of the two types of 
isotropy subgroups as 

H^il^o)) = "i?7 ° ^^IV'o) " = {76 Tp s.t. pexp(7) = h, with h G (13) 

and call H^{\ipo)) the ^-isotropy subgroup at \ipo)- 

In system theory terminology, H^Qipo)) should be called the "unobservable subspace" of the "state space", 
because "states"_7 G H^{\iIjq)) are indistinguishable when looked from the "output" \tp)- 

Similarly to (y), one can define the -f-coset space Cj{\4>q)) as 

C7(IV'o)) = rp/i/,(|^o)) 

and call the ^-coherent states the values assumed by \ip) for 7 varying in C-y ( | ■00 ) ) ■ 
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B. Universality of the gates for products of exponentials coordinates 

The universality of the gates property mentioned in the introduction corresponds to C^dipo)) being all of 
giv-i^ holds almost everywhere and depends on the free Hamiltonian (i.e. on the parameters a.;) being 
nondegenerate, see [Q. Notice that in the present framework the correct system theoretic concept to use is 
"output" controllability, not "state" controllability. Notice, furthermore, that whenever S C H^{\iPq)), the 
singular points of the "state" space representation (|l0| ) are ininfluent in the quantum state manipulation, 
i.e. they are not seen from the "output" 1-0) of (|lO|). Thus we have the following result: 

Proposition 1 Assume that the system (|^) is controllable and that a fixed order of infinitesimal generators 
is chosen in (^). For a quantum state \ipo), 

i) if ^ Q Hydtpo)), then the gates in (|^) are universal; 

ii) if Ti H~^(\^q)), then the gates in ^ may not be universal. In this case the analysis of the singular 
locus Yi is required in order to establish universality of the gates. 

While the value of S depends only on the order chosen in (^, iJ^(|-0o)) varies also with the initial condition 
I'i/'o). Thus the previous sufficient condition for universality of the gates may hold only for some values 
of I'i/'o)- Notice that, provided that (|l|) is controllable, the sufficient conditions above does not depend on 
the control authority available (i.e. on how many Uj, j = 1, . . . , n are different from 0). However, when 
S ^ Hydipo)) and the singular points of E are investigated, universality depends essentially on which control 
fields are available. 

All the constructions and characterizations discussed so far are easily clarified by means of a concrete 
example. 



IV. APPLICATION TO A SINGLE QBIT 



The most elementary quantum state of interest here is the single qbit, defined on the sphere S^. In the 



computational basis |0}, |1}, the state is — [yi J/2] with 
is expressed in the computational basis coordinates as 



\y2\ 



1, thus the complex sphere in 



§1 



{ 


yi 




V2 



SU{2) acts transitively on and we shall concentrate on SU{2) unitary transformations as elementary gates 
for \ip). A skew-symmetric basis for su(2) is obtained from the Pauli matrices 





1 




-i 




'1 " 




1 


(^2 = 


i 


<^3 = 


-1 



for example by taking Aj — -^(Jj, j = 1, 2, 3, i.e. 



Ai = - 
2 



i 

1 



1 

-1 



A. = ^ 



i 
-i 



(14) 



and it corresponds to all real structure constants — — Cgj^ = 1. The corresponding adjoint matrices 
are 



ad^i 

whose exponentials are already known from the literature to be: 



"0 





" 




" 





1" 




"0 


-1 0' 








-1 


ad^2 = 











adA3 = 


1 








1 







-1 

















1 

COS 71 — sin 71 
sin 71 COS 71 



COS 72 sin 72 

1 
- sin 72 COS 72 



COS 73 — sin 73 
sin 73 COS 73 
1 
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The parameter space T i.e. the portion of the Lie algebra that is needed to parameterize the group under 
the single exponential is a solid sphere of radius 27t and with all the points on the surface identified (they 
all correspond to —I, see 0, p. 123). 



A. ZYZ quantum logic gates 

Probably the most common set of gates for the single qbit is given by the ZYZ Euler angles. Choosing 
such a parameterization for SU{2), the ZYZ-operations on lijjo) are described by the following product of 
exponentials 



\^) = (7(71, 72 , 73)|V'o) = e'^i^e'^-^e'^^^li^o) = e'<'^'e'<'^- e-'-^'lyjo) 



(15) 



where 7 = [71 72 73] is defined in Fp ~ (— Stt, 27r]'^ and 7^ = %{t). 
From the expressions for the exponentials of the Pauli matrices: 



o73^3 _ 



7l r 

COS-/ 

COS-/ 

73 r 
COS — / 



. . 7i 
I sm "^o'l 

. . 72 

I sm — CT2 
. . 73 

I sm — CTs 
2 



cos ^ i sin ^ 

21 r-na 21 



cos sm 



e'T- 

- "7 

e-*^ 



the ZYZ product of exponentials is 



U{-fi, 72, 73) 



cos ■ 



2 ' sm ■ 



2 sm - 

23-1 

2 COS 



(16) 



B. Wei- Norman formula for the ZYZ Euler angles 



The Wei-Norman formula corresponding to the ZYZ ordered product (|T^) is (see 



and its inverse 



—sin 71 cos 71 sin 72 

cos 71 sin 71 sin 72 

1 cos 72 



— cos 7i cot 72 — sin 71 cot 72 1 
— sin 71 cos 71 

cos 71 CSC 72 sin 71 esc 72 



(17) 



(18) 



Since 



det(S) — sin 72 

fcTT, fc G Z, as is well-known for such a parameterization. Thus, in Fp 



the singular points correspond to 72 
E!~^ can be used everywhere except in S = {7 S Fp s.t. 72 = — tt, 0, tt, 27r}. 

It is worth emphasizing that it is a fundamental topological fact that singularities cannot be avoided 
in a minimal parameterization of a semisimple Lie group. One possible way to get around the problem is 
obviously to use "redundant" parameterizations like quaternions, another to choose two different local charts 
i.e. two sets of of exponential coordinates corresponding for example to a different ordering in (15). 
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C. Example: nuclear spin qbit 

In (pO|), some of the aj and/or Uj might be zero. If for example the free Hamiltonian Hq is diagonal then 
fli = a2 = 0. In this case, most likely U3 = if the dipole approximation for the coupling with the control 
fields is considered. Assume that independent fields along both X and Y directions are available. Using (p^), 
(p^ reads as 



71 




as 




— cos 71 cot 72 


— sin 71 cot 72 


72 







+ 


— sin 7i 


COS 7i 


73 









cos 71 CSC 72 


sin 71 CSC 72 



Vi 




V2 





' COS ^ 



o/v 2^2) sin — 



) sin ^ 
COS ■ 



(19) 



2/I0 



The model ( |l9|) corresponds to the well-known example of a single nuclear spin ^ system immersed in a 
magnetic field which is constant along the Z axis and varying in the X and Y directions. 



D. Coherent states for 5(7(2) 
The isotropy subgroup of a generic state \4>) can be obtained by solving the algebraic equation 

|V') = C/(7i,72 , 73)|V') (20) 

The ZYZ representation is convenient is this respect since it is easy to notice that (^0|) is solved for all \ipo) at 
least by the submanifold 72 = and 73 = —71. In fact C/(7i, 0, —71) = / and therefore Hj{\^po)) ^ H^{I) = 
{7 G Fp s.t. 72 = 0, 73 = —71}. A complete computation of the Hjdtpo)) depends on the value of lipo)- For 
example, for |0) ~ [1 0]^ ( ^ ) is solved by all 7 such that 72 = or 72 = 27r, since in this case e^^^'*'"^'^ 
becomes an ininfluent global phase factor: |0) = ±e*(^+^)|0). Thus H^{\0)) {7 g Fp s.t. 72 = 0, 27r}. 
For 72 = and ^ + ^ = f , this corresponds to the so-called "Z gate": C/(7i, 0, | -71) = e'^i^^e^^-'^^^^^ = 



1 

-1 

l^. If instead {ipo) is a superposition {yi^ ^ and j/2o ^ 0), then the 7-isotropy subgroup contains H^{I) 
but may also contain other points, according to the value of IV'o)- We are in the case ii) of Propositior 
S ^ iJ^d^/^o)) for some jV'o) (other than |0) and |1)). 



In fact |0) is a highest weight state in for which the computation of coherent states is easier 

^(7) 



E. Analysis of the singularity of the Wei-Norman formula 

In E, the Wei-Norman formula (|l8|) is not valid anymore. However, in this case (^^ becomes: 



5(71, 0, 73) 



— sin 71 

cos 71 

1 1 



(21) 



Since in S, rank 2(71, 72, 73) = 2, varying u at most two of the 7^ can be varied independently. From ( pl| ) 
and (0) we get 



71+73 = 03 + 1*3 

72 = — (ai + ui) sin7i + (a2 + -^2) COS71 



(22) 
(23) 



i.e. it is still possible to obtain a well-defined expression for the dynamics at which 72 obeys. 
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In particular, in the example of (|l9|), the system in S reduces to 



71 + 73 

72 



as 

-Ui sin 7i + U2 cos 71 



-/ 11+13 \ 

_,Y21±23 \ 

e ^ 2 •'y2o 



from which it is easy to observe that the free Hamiltonian only induces a relative phase transformation on a 
superposition state \4>o) (as expected, Hq being diagonal), while the actuated system can always be steered 
out of the "unobservable subspace" H^{\tpo)) (and also of E) by means of the available control authority ui 
and U2- 

In summary: almost all the singular points of S are "unobservable" for \ip), i.e. they produce no influence 
on the quantum state. Some care is required in the algebraic set S because the exponential coordinates are 
not defined. The problem is analogous to the lack of global coordinates on a manifold. Even if is not 
defined in E, the implicit equation (||) still governs the dynamics and can be used to pull 7 out of E. 



V. CONCLUSION AND OUTLOOK 



The aim of this work is to build the differential equations of a driven quantum dynamics in the formalism 
best suited for quantum computation, where a predefined set of quantum gates can be expressed directly 
in terms of the physical control parameters via exponential coordinates on the group of transformation. 
Beside the general remark that choosing coordinates on a manifold is always necessary when one wants to 
do computations, in the present context it is worth stressing that coordinates (here exponential coordinates) 
have the appealing interpretation of "gains" in the logic gates of the quantum hardware. 

We would like to emphasize that the method presented here is by no means the only way to relate the 
continuous and discrete dynamics. Indeed each implementation proposal for quantum computers has its own 
method to generate gates, relying on the already existing techniques for state manipulation. See fl^ for a 
survey of the different methods. The practice on NMR spectroscopy for example, bypasses the synthesis of 
the control fields in continuous time by applying strong pulses of predetermined time width and by checking 
the overall effect on the Hamiltonian via coherent averaging and syncronous stroboscopic measurements Q . 
The use of piecewise constant controls mentioned in the Introduction is indeed another possible solution. 
On the contrary, notice that the well-established formalism of sampled (bilinear) control systems, which 
one could think of using in the present context, is actually not suited, as the potential power of quantum 
computing resides in the continuous superposition of quantum states which is destroyed by sampling. 

The method we are proposing here not only enable to analytically describe the unitary transformations 
in terms of the available control inputs, but it also gives the possibility of monitoring their variations due 
to dynamical perturbations or uncertainties in the external fields. Since we choose to work with exponential 
coordinates for the unitary group SU(N) of the TVqbit, non unitary disturbances like decoherence phenomena 
cannot be represented. We plan to extend the method to open quantum systems in the future. 
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